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Abstract 

The dynamical mass generation for gluons is discussed in Euclidean Yang-Mills theories 
supplemented with a renormalizable mass term. The mass parameter is not free, being 
determined in a self-consistent way through a gap equation which obeys the renormalization 
group. The example of the Landau gauge is worked out explicitly at one loop order. A few 
remarks on the issue of the unitarity are provided. 
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1 Introduction 



In the last years many efforts have been done to put forward the idea that gluons might ac- 
quire a mass through a dynamical mechanism. These efforts have led to a considerable amount 
of evidence, obtained through theoretical and phenomenological studies [H EJ El HJ EJ EJ El 
El El EB HIl OS US1 IIZl QE1 QS1 HQ1 as well as from lattice simulations 

[231 EH EH EH EH EH1 E9l E01 ED E21 [33]. Many aspects related to the dynamical gluon mass 
generation deserve a better understanding. This is the case, for example, of the unitarity of the 
resulting theory, a highly nontrivial topic, due to the confining character of QCD. 

Needless to say, the unitarity of the S matrix is a fundamental property of the spectrum of 
a quantum field theory. It expresses the conservation of the probability of the amplitudes cor- 
responding to the various scattering processes among the excitations of the spectrum. 

In a nonconfining theory, the first step in the construction of the S matrix is the introduc- 
tion of the so-called \in) and \out) Fock spaces characterizing the asymptotic behavior of the 
physical states before, t — ► — oo, and after, t — > +00, a scattering process. The S'-matrix is thus 
defined as the unitary operator which interpolates between the spaces \in) and \out), namely 

I in) = S \out) . (1) 

The relation of this equation with the Green's functions of the theory is provided by the LSZ 
formalism. A key ingredient of this formalism is the introduction of the asymptotic fields, ipi n 
and (fout, describing the asymptotic behavior of the interacting fields ip, according to 

¥>lt— 00 = ^W, (2) 
Plt-H-oc = Z "Pout ■ 

The asymptotic fields cpi n and <p ut allow us to define the creation and annihilation operators 
( a L> a *«) an d ( a Ltt> a autj, from which the Fock spaces \in) and \out) are obtained. The entire 
construction relies on the possibility that the asymptotic fields can be consistently introduced. 



In a confining theory like QCD, the quanta associated with the basic fields of the theory, 
i.e. the gluon field and the quark fields ifi, ip, cannot be observed as free particles, due to 
color confinement. The physical spectrum of the theory is made up by colorless bound states 
of quarks and gluons giving rise, for instance, to barions, mesons and glueballs. This implies 
that the asymptotic Fock spaces \in), \out) of the theory have to be defined through suitable 
operators from which the physical spectrum of the excitations is constructed. Of course, the 
S'-matrix describing the scattering amplitudes among the excitations of the physical spectrum 
of the theory has to be unitary. Although intuitively simple and easily understandable, this 
framework is far beyond our present capabilities. An operational definition of the gauge invari- 
ant colorless operators defining the physical spectrum of the excitations and a well defined set 
of rules to evaluate their scattering amplitudes are not yet at our disposal. 

Quantized Yang-Mills theories are described by the Faddeev-Popov LagrangiarQ 

S = S YM + S gf = J d A x (jF%F%, + b a d^ + c% (D^cA , (3) 

here taken in the Landau gauge. The field b a in expression ^) is the Lagrange multiplier 
enforcing the Landau gauge condition, d^A^ = 0, while c a , c a stand for the Faddeev-Popov 

*From now on we shall consider pure Yang-Mills theory in the Euclidean space-time. 



2 



ghosts. The action (|3|) is renormalizable to all orders of perturbation theory and displays color 
confinement^. Furthermore, thanks to the asymptotic freedom, the gauge field behaves 
almost freely at very high energies, where perturbation theory is reliable. However, at low 
energies, the coupling constant increases and the effects of color confinement cannot be neglected. 
We do have thus a good understanding of the properties of the field A® at high energies, whereas 
it becomes more and more difficult to have a clear picture of A^, and of the whole theory, as the 
energy decreases. We might thus adopt the point of view of starting with a renormalizable action 
built up with a gauge field A^ which accommodates the largest possible number of degrees of 
freedom. This would amount to start with a quantized massive Yang-Mills action 

S m = SyM + Sgf + Smass , (4) 

where S mass is a suitable mass term for the gauge field A^. The best choice for S mass would 
be a gauge invariant, renormalizable local mass term. However, no local renormalizable gauge 
invariant mass term built up with gauge fields only is at our disposal. Nevertheless, it might be 
worth reminding that, recently, a consistent framework for the nonlocal gauge invariant mass 
operator 
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ab 



0(A) = -- I #xF%, \P ) F%. (5) 



has been achieved [31] . More precisely, the nonlocal operator (0) can be cast in local form 
by means of the introduction of a suitable set of additional fields. The resulting local theory 
displays the important property of being multiplicatively renormalizable 



Though, for the time being, we give up of the requirement of the gauge invariance. This will 
enable us to present our analysis with the help of a relatively simple example. Therefore, as 
possible mass term we shall take 

Smass = \m 2 I #xA%A% , (6) 
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so that 



S m = J d 4 x {jF; v F; v + l m 2 AlAl + b a d^Al + lfd,{D^ . 



(7) 



Expression ([7]) provides an example of a massive nonabelian gauge theory which is renormaliz- 
able to all orders of perturbation theory [35] , while obeying the renormalization group equations 
(RGE). 

A few remarks are now in order: 



• The amplitudes corresponding to the scattering processes among gluons and quarks display 
now a violation of the unitarity. This can be understood by noting that the inclusion of 
the mass term mP'A'tA't gives rise to a BRST operator which is not nilpotent. However, 
as shown in [35j, it is still possible to write down suitable Slavnov- Taylor identities which 
ensure that the massive theory (|7|) is renormalizable to all orders of perturbation theory. 
Moreover, if sufficiently small, this violation of the unitarity might not be in conflict 
with the confining character of the theory. Otherwise said, since gluons are not directly 
observable, we could allow for a gauge field A a ^ with the largest possible number of degrees 
of freedom, provided the renormalizability is preserved and one is able to recover the results 
of the massless case at very high energies. 

* Although we still lack a definite proof of color confinement, it is doubtless that pure Yang-Mills theory, as 
given in eq.Q, displays such a phenomenon. 
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• This framework would be useless if the value of the mass parameter m would be free, 
meaning that we are introducing a new arbitrary parameter in the theory, thereby chang- 
ing its physical meaning. A different situation is attained by demanding that the mass 
parameter is determined in a self-consistent way as a function of the coupling constant g. 
This can be obtained by requiring that the mass m in eq.(J7|) is a solution of a suitable gap 
equation. In other words, even if the mass m is included in the starting gauge-fixed theory, 
it does not play the role of a free parameter, as it is determined once the quantum effects 
are properly taken into account. Here, we rely on the lack of an exact description of a con- 
fining Yang-Mills theory at low energies. We start then with the largest possible number 
of degrees of freedom compatible with the renormalizability requirement and fix the mass 
parameter through the gap equation. If the resulting value of m will be small enough, one 
can argue that the unitarity is violated by terms which become less and less important 
as the energy of the process increases, so that the amplitudes of the massless case are in 
practice recovered at very high energies. The present set up might thus provide a different 
characterization of the aforementioned phenomenon of the dynamical gluon mass genera- 
tion, which has already been successfully described in [U [TJ [TUl ttH [J3J [TU [IBJ E]- ^ n ^ ne 
next section, the gap equation for the mass m will be discussed. 



2 The gap equation for the mass parameter m 

The gap equation for the mass parameter m is obtained by requiring that the vacuum functional 
£ defined by 

~ V£ = I [23*] e ~ t) , (8) 



where V is the Euclidean space-time volume, obeys a minimization condition with respect to 
the mass m, i.e. the value of the mass m is determined by demanding that it corresponds to 
the minimum of the vacuum functional £, namely 

8£ 0. (9) 



dm 2 



Equation ^ is the gap equation for the mass parameter m. The quantity rj(g) in eq.([S]) is a 
dimensionless parameter whose loop expansion 

v(g) = m(g) + + h 2 m{g) + (io) 

accounts for the quantum effects related to the renormalization of the vacuum diagrams in the 
massive case. The parameter r](g) can be obtained order by order by requiring that the vacuum 
functional £ obeys the renormalization group equations (RGE) 

d£ , . 

"rr"' (11> 

meaning that £ is independent from the renormalization scale /U, as it will be explicitly verified 
in the next section. Equation (jllf) expresses an important property of the vacuum functional 
£. We also remark that a term of the kind of 77m 4 in eq.([8l) has been already obtained^ in [1] in 
the evaluation of the vacuum energy of Yang- Mills theories when gluons are massive. 

The gap equation equation (llip can be given a simple interpretation. Due to the lack of an 



*See eq.(6.22) of Sect.VI of [I]. 
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exact description of Yang-Mills theories at low energies, we have adopted the point of view of 
starting with a renormalizable massive action, as given in eq.©. As far as the mass parameter to 
is free, expression ([7]) can be interpreted as describing a family of massive models, parametrized 
by to. For each value of m we have a specific renormalizable model. Moreover, as the introduc- 
tion of a mass term has an energetic coast, we might figure out that, somehow, the dynamics 
will select precisely that model corresponding to the lowest energetic coast, as expressed by the 
gap equation (fTTI) , 

Before starting with explicit calculations let us summarize our point of view: 



• Since gluons are not directly observable, we allow for a gauge field with the largest 
number of degrees of freedom compatible with the requirement of renormalizability. 

• This amounts to start with a renormalizable massive action, as given in eq.([7]). However, 
the mass parameter to is determined in a self-consistent way by imposing the minimizing 
condition ([9]) on the vacuum functional £. 

• Also, it is worth observing that, in the case of the massive model of eq.(J7|), a non vanishing 
solution, m? sol ^ 0, of the gap equation ([9]) implies the existence of a non vanishing dimen- 
sion two gluon condensate (A®Afy. In fact, differentiating equation ([8]) with respect to 
to 2 and setting to 2 = m 2 ol , one obtains 

l(A a ^) = - V m 2 sol . (12) 
3 Evaluation of the vacuum functional 8 at one loop order 

In the case of pure SU(N) Yang-Mills theories, for the vacuum functional £ we have 

e -™ = J m A Sm+Vr >^) , (13) 
with S m given by expression (J7j), namely 

S m = J d±x (±F%F% + \m % AlAl + b a d,A« + c a d, (D^cA . (14) 



As it has been proven in [35J, the massive action (114j) is multiplicatively renormalizable to all 
orders of perturbation theory. In particular, for the mass renormalization we have [35] 

90 = Z g g , 
A = Z X l 2 A 

2 ry 2 

TOq = Z m 2TO , 

Z m 2 = Z g Z A I , (15) 



from which the running of the mass to 2 is easily deduced 

dm- 
d/j, 



dm 2 2 

M~5 = -7m2TO, , (16) 



with 

i m >(g 2 ) = iog 2 + iig A + 0{g 6 ) , (17) 
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35 N 449 / N x 2 



70 = yi6^' 7i = ^U-v • (18) 

Also 

P(g 2 ) =JI^ = -2 (A,/ + Pig 6 + 0(g 8 )) , (19) 
11 N „ 34 / N x 2 



*-sW A = Tll6^J ' (20> 

In order to obtain the parameter rj at one-loop order, it is useful to note that expression (|13p 
can be written in localized form as 



e 



vt - J e -( 5 ™ +y ^) = | DJ(x) 5(J(x) - m 2 ) e~ w ^ , (21) 



with 

-W{J) 



e 



J [Z>$] e"^ , (22) 
S-(J) = 5yM + < S 9/ + y' ( i 4 xQj( 2 ;)^ + |j 2 (x)) . 



From equation (|2ip it follows that the renormalization of the vacuum functional 8 can be 
achieved by renormalizing the functional W(J) in the presence of the local source J(x), and 
then set J = m? at the end. The renormalization of the functional W(J) has been worked out 
at two-loops in [7J. By simple inspection, it turns out that the parameter rj is related to the 
LCO parameter £ of [7] by rj = yielding 

9 N 2 — 1 161 iV 2 - 1 2 



Thus, for the vacuum functional £ at one-loop order in the MS scheme, we get 

„ m 4 / 9 iV 2 -l 161iV 2 -l\ ^iV 2 -l A ( 5 , m 2 \ . jN 

where we have introduced the factor h to make clear the order of the various terms. It is useful 
to check explicitly that the above expression obeys the RGE equations. Indeed, from eqs. (!17p . 
(fTUj) we obtain 

- d£ * 2 if 9 N 2 -l\ ,m A 9 iV 2 -l _ 4 _iV 2 -l 4 , 

= n m ^(^)^-n m ^ 3 l-m^l^ + (n 2 ) 

16vr 2 V 6 / 13 167r 13 647r 

fc 4 iV 2 -l/35 9 33 6\ „,_ 2 . . 4 N 2 - 1 /105 33 3\ _ 2 , 



16tt 2 V 6 13 13 A J K ' 16vr 2 V 26 I 3 2 # 

, 4 iV 2 -l /105 - 66 - 39\ 2 2 

= ^Sfo^V 26 J+0(^) = 0(^). (25) 

It remains now to look for a sensible solution of the gap equation © . This will be the task of 
the next section. 
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3.1 Searching for a sensible minimum 



In order to search for a sensible solution of the gap equation Q, -j^p = 0, we first remove the 
freedom existing in the renormalization of the mass parameter by replacing it with a renormal- 
ization scheme and scale independent quantity. This can be achieved along the lines outlined in 
|10| in the analysis of the gluon condensate (A®Afy within the 2PPI expansion technique. Let 
us first change notation 

g 2 -> f , (26) 

o 2 

m — > m , 



and rewrite the one-loop vacuum functional as 



S 



9 iV 2 - 1 1 
13 N ]? 



m 4 13 Ng 2 ,/, m 2 113 
m lor 



(27) 



2 ' 3 64vr 2 V J? 39 
As done in [10], we introduce the scheme and scale independent quantity m 2 through the relation 

m 2 = J(g 2 )m 2 . (28) 

From 

77- 



¥^ = -lm<t)™\ (29) 



with 

lm<t) = 709 2 + 7i5 4 + O(5 6 ) , (30) 

35 N _ 449 f N \ 2 

70 = yi6^' 7i = ^vi6^J ' (31) 



we obtain the condition 



from which it follows that 



H^ 1 =^(9 2 mf), (32) 



djl 



dp 

Equation (|32|) is easily solved, yielding 



^ = . (33) 



f(f) = (g 2 T^{i + hg 2 + o(t)) , 

fo = ^r(^i-ii) , (34) 



2po V/3o 



where the coefficients /3q, 0i are given in eqs. (|19p . (|2Q|) . Moreover, one has to take into account 
that a change of scheme entails a change in the coupling constant g 2 , according to 

g 2 = g 2 (l + b g 2 + O(g 4 )). (35) 

The coefficient bo in eq. (|35p expresses the freedom related to the choice of the renormalization 
scheme. It will be fixed by demanding that the coupling constant is renormalized in such a 
scheme so that the vacuum functional 8 takes the form 

„,_ 2 . 9iV 2 -l 1 / m 4 ^ 4 Ng 2 ^ \ 
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where L stands for 



L = log 



m 



2 (g 2 )^ 



~p 2 



(37) 



and Ei is a numerical coefficient. After a simple calculation, we get 



9 N 2 



r ~ 4 
m 



13 N 



Or 



^„ 4 13Ng 2 { 113 3 64vr 2 / 6q n 7(K 



Therefore, for 6q ° ne nas 



namely 



113 3 647T 2 / 6o n 7ck 



39 



0, 



bo 



4331 TV 
" 396 16vr 2 



For the vacuum functional £ (m 2 ) one gets 

„ 9 iV 2 - 1 1 



13 N 



(9 2 t 



70 



m 4 ^ 4 13A^ 2 

~T + m Y64^ 



(38) 

(39) 
(40) 

(41) 



In terms of the scale independent variable m 2 , the gap equation reads 

d£ 



dm 2 



0, 



so that 



~ 2 ~ 2 

— m + m 



26/Y 5 2 2 13iV 5 2 
r^riy + m — T = . 



3 64vr 2 



3 64vr 2 



Next to the solution, m 2 = 0, we have the nontrivial solution m so i given by 



7Q 



3 64vr 2 



7J 2 



13 Ng 2 n 
-I — = o 

^ 3 64vr 2 



(42) 
(43) 

(44) 



In order to find a sensible solution of this equation, a suitable choice of the scale fx has to be 
done. Here, we take full advantage of the RGE invariance of the vacuum functional £ , and set 



f 1 =m S oi 



(g 2 ) 2/3 ° 



(45) 



where s is an arbitrary parameter which will be chosen at our best convenience. The possibility 
of introducing this parameter relies on the independence of the vacuum functional £ from the 
renormalization scale Ji. Furthermore, recalling that 



g 2 m 



(46) 



and that, due to the change of the renormalization scheme, 



A 



A- e 0o 



(47) 
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for the effective coupling and the mass m so i, one finds 



Ng 2 



12 



1 



16vr 2 
12 16vr 2 



l— loop 

1 



13 (1 + 2s) 



7Q 
"4/3 



e 2/3o e S( 1 + 2s )ef A 7 



m ^ll-ioop- \^ 13 JV (1 + 2s) — ■ — M5- 

Therefore, choosing s = 0.6, and setting N = 3, the following one-loop estimates are found 

2 

1— Zoop 



16vr 2 



0.42 , 



m sol\i-l oop 



2 4A ~ 

233iWW , 



560MeV 



(48) 
(49) 

(50) 
(51) 



7V=3 



1— Zoop 



0.22GeF 



and 



c /~ \iJV=3 

^("Wjli-joop 



-90A- ~ 



-0.265 (GeF) 4 



(52) 



Note that the value obtained for fh so i is close to that already reported for the dynamical gluon 
mass in the Landau gauge [TOj [HI [13 E3 EI]- It should be remarked that the results ([50]) . 
(|5ip have been obtained within a one-loop approximation. As such, they can be taken only 
as a preliminary indication. To find more reliable results, one has to go beyond the one-loop 
approximation. Nevertheless, these calculations suggest that a non vanishing gluon mass might 
emerge from the gap equation Q. 



4 Conclusion 

In this work the issue of the dynamical mass generation for gluons has been addressed. Due to 
color confinement, gluons are not observed as free particles. Thanks to the asymptotic freedom, 
the gauge field behaves almost freely at very high energies, where we have a good under- 
standing of its properties. However, as the energy decreases the effects of confinement cannot 
be neglected and it becomes more and more difficult to have a clear understanding of A^. As 
a consequence, one does not exactly know what is the correct starting point in the low energy 
region. We might thus adopt the point of view of starting with a renormalizable action built up 
with a gauge field A^ which accommodates the largest possible number of degrees of freedom. 
This would amount to take as starting point a renormalizable massive action, as considered, for 
example, in expression (J3). The mass parameter m is not treated as a free parameter. Instead 
it is determined by a gap equation, eq.(|9|), obtained by minimizing the vacuum functional 8 of 
eq.(Hl) with respect to the mass parameter m. A preliminary analysis of this gap equation at 
one-loop shows that a nonvanishing gluon mass might emerge. Also, the vacuum functional £ 
displays the important feature of obeying the renormalization group equations. 

Finally, we underline that the infrared behavior of the gluon propagator is expected to be 
affected by several mass parameters, with different origins. For instance, as pointed out in 
[361 137] in the case of the Landau gauge, the gluon propagator turns out to be affected by both 
dynamical gluon mass m and Gribov parameter 7, which arises from the restriction of the do- 
main of integration in the Feynman path integral up to the first Gribov horizon. More precisely, 
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these parameters give rise to a three level gluon propagator which exhibits infrared suppression 
157] , namely 



Acknowledgments. 



I am indebted to my friends and colleagues D. Dudal, J. A. Gracey, and H. Verschelde for 
many valuable discussions. T. Turner and L. R. de Freitas are gratefully acknowledged. We 
thank the Conselho Nacional de Desenvolvimento Cientffico e Tecnologico (CNPq-Brazil), the 
Faperj, Fundagao de Amparo a Pesquisa do Estado do Rio de Janeiro, the SR2-UERJ and the 
Coordenagao de Aperfeicoamento de Pessoal de Nivel Superior (CAPES) for financial support. 



References. 

[1] J. M. Cornwall, Phys. Rev. D 26, 1453 (1982). 

[2] J. Greensite and M. B. Halpern, Nucl. Phys. B 271, 379 (1986). 

[3] M. Stingl, Phys. Rev. D 34, 3863 (1986) [Erratum-ibid. D 36, 651 (1987)]. 

[4] M. J. Lavelle and M. Schaden, Phys. Lett. B 208, 297 (1988). 

[5] F. V. Gubarev and V. I. Zakharov, Phys. Lett. B 501, 28 (2001) jarXiv:hep-ph/0010096] . 

[6] F. V. Gubarev, L. St odolsky and V. I. Zakharov, Phys. Rev. Lett. 86, 2220 (2001) 
[arXiv:hep-ph/0010d57| . 

[7] H. Verschelde, K. Knecht, K. Van Acoleyen and M. Vanderkelen, Phys. Lett. B 516, 307 
(2001) | arXiv:hep-th/0105018| . 

[8] K. I. Kondo, Phys. Lett. B 514, 335 (2001) [arXiv:hep-th/0105299] . 

[9] K. I. Kondo, T. Murakami, T. Shinohara and T. Imai, Phys. Rev. D 65, 085034 (2002) 
[arXiv:hep-th/01112~56] . 

[10] D. Dudal, H. Verscheld e, R. E. Browne and J. A. Gracey, Phys. Lett. B 562, 87 (2003) 
[arXiv:hep-th/0302128] . 

[11] R. E. Browne and J. A. Gracey, JHEP 0311, 029 (2003) [arXiv:hep-th/0306200| . 

[12] U. Ellwanger and N. Wschebor, Int. J. Mod. Phys. A 18, 1595 (2003) 
[arXiv:hep-th/0205057] . 

[13] D. Dudal, H. Verschelde, V. E. R. Lemes, M. S. Sarandy, S. P. Sorella and M. Picariello, 
Annals Phys. 308, 62 (2003) |arXiv:hep-th/0302168] . 

[14] D. Dudal, H. Verschelde, J. A. Gracey, V. E. R. Lemes, M. S. Sarandy, R. F. Sobreiro and 
S. P. Sorella, JHEP 0401, 044 (2004) [arXiv:hep-th/0311194| . 

[15] A. C. Aguilar, A. A. Natale and P. S. Rodrigues da Silva, Phys. Rev. Lett. 90, 152001 
(2003) [arXiv:hep-ph/0212105] . 



10 



[16] D. Dudal, J. A. Gracey, V. E. R. Lemes, M. S. Sarandy, R. F. Sobreiro, S. P. Sorella and 
H. Verschelde, Phys. Rev. D 70, 114038 (2004) [arXiv:hep-th/0406132| . 

[17] R. E. Browne and J. A. Gracey, Phys. Lett. B 597, 368 (2004) [arXiv:hep-ph/0407238] . 

[18] J. A. Gracey, Eur. Phys. J. C 39, 61 (2005) jarXiv:hep-ph/0411169] . 

[19] G. Parisi and R. Petronzio, Phys. Lett. B 94, 51 (1980). 

[20] J. H. Field, Phys. Rev. D 66, 013013 (2002) [arXiv:hep-ph/0101158| . 

[21] A. P. Szczepaniak and E. S. Swanson, Phys. Lett. B 577, 61 (2003) [arXiv:hep-ph /0308268] . 

[22] X. d. Li and C. M. Shakin, Phys. Rev. D 71, 074007 (2005) [arXiv:hep-ph/0410404] . 

[23] K. Amemiya and H. Suganuma, Phys. Rev. D 60, 114509 (1999) [arXiv:hep-lat/9811035| . 

[24] P. Boucaud, A. Le Yaouanc, J. P. Leroy, J. Micheli, O. Pene and J. Rodriguez-Quintero, 
Phys. Rev. D 63, 114003 (2001) [arXiv:hep-ph/0101302| . 

[25] P. Boucaud et al, Phys. Rev. D 66, 034504 (2002) [arXiv:hep-ph/0203119] . 

[26] Ph. Boucaud et al, Phys. Rev. D 74, 034505 (2006) [arXiv:hep-lat /0504017| . 

[27] K. Langfeld, H. Re inhardt and J. Gattnar, Nucl. Phys. B 621, 131 (2002) 
[arXiv:hep-ph/010714l] . 

[28] C. Alexandrou, P. de Forcrand and E. Follana, Phys. Rev. D 63, 094504 (2001) 
[arXiv:hep-lat/00080"l2] . 

[29] C. Alexandrou, P. de Forcrand and E. Follana, Phys. Rev. D 65, 114508 (2002) 
[arXiv:hep-lat/01 12043] . 

[30] V. G. Bornyakov, M. N. Chernodub, F. V. Gubarev, S. M. Morozov and M. I. Polikarpov, 
Phys. Lett. B 559, 214 (2003) [arXiv:hep-lat/0302002] . 

[31] E. Ruiz Arriola, P. O. Bowman and W. Broniowski, Phys. Rev. D 70, 097505 (2004) 
[arXiv:hep-ph/0408309l . 

[32] T. Suzuki, K. Ishiguro , Y. Mori and T. Sekido, Phys. Rev. Lett. 94, 132001 (2005) 
[arXiv:hep-lat/0410001~l . 

[33] M. N. Chernodub et al, Phys. Rev. D 72, 074505 (2005) [arXiv:hep-lat /0508004| . 

[34] M. A. L. Capri, D. Dudal, J. A. Gracey, V. E. R. Lemes, R. F. Sobreiro, S. P. Sorella and 
H. Verschelde, Phys. Rev. D 72, 105016 (2005) l iarXiv:hep-th/0510240| . 

[35] D. Dudal, H. Versc helde and S. P. Sorella, Phys. Lett. B 555, 126 (2003) 
[arXiv:hep-th/0212182] . 

[36] R. F. Sobreiro, S. P. S orella, D. Dudal and H. Verschelde, Phys. Lett. B 590, 265 (2004) 
[arXiv:hep-th/0403135] . 

[37] D. Dudal, R. F. Sobreir o, S. P. Sorella and H. Verschelde, Phys. Rev. D 72, 014016 (2005) 
[arXiv:hep-th/0502183] . 



11 



